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Abstract 

With the aim of constraining the structural properties of neutron stars and the 
equation of state of dense matter, we study sudden spin-ups, glitches, occurring 
in the Vela pulsar and in six other pulsars. We present evidence that glitches 
represent a self-regulating instability for which the star prepares over a waiting 
time. The angular momentum requirements of glitches in Vela indicate that 
> 1.4% of the star's moment of inertia drives these events. If glitches originate 
in the liquid of the inner crust. Vela's 'radiation radius' i?oo must exceed ~ 
12 km for a mass of IAMq. The isolated neutron star RX J18563-3754 is a 
promising candidate for a definitive radius measurement, and offers to further 
our understanding of dense matter and the origin of glitches. 



INTRODUCTION 

Many isolated pulsars suffer spin jumps, glitches, superimposed upon oth- 
erwise gradual spin down under magnetic torque. For example, in the glitch of 
the Crab pulsar shown in Fig. 1, the star spun up by nearly a part in 10'' over 
several days [1], corresponding to a change in rotational energy of the crust of 
~ 10^^ ergs. A particularly active glitching pulsar is the Vela pulsar, which has 
produced more than a dozen glitches since its discovery over 30 years ago. The 
fractional changes in rotation rate are typically ~ 10~^, occurring every three 
years on average [2]. 

Because glitching pulsars are isolated systems, glitches are thought to arise 
from internal torques exerted by the rotating liquid interior on the crust, whose 
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spin rate we observe [3]. As the star's crust is spun down by the magnetic field 
frozen to it, the interior liquid, which responds to the external torque indirectly 
through friction with the solid crust, rotates more rapidly. For example, a portion 
of the liquid could coast between glitches while the solid crust spins down (Fig. 
2). Glitches might arise as the consequence of an instabiUty that increases the 
frictional coupling between the liquid and the solid, causing angular momentum 
flow to the crust. 
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Figure 1 A glitch in the Crab pulsar. The 
spin rate with respect to a pre-glitch spin- 
down model is shown. Data are from Ref. 
1. 



Figure 2 A glitch model. Differential rota- 
tion between the interior liquid (dashed line) 
and the crust (solid line) develops as the crust 
is spun down by its external torque. At the 
time of a glitch, the liquid delivers angular 
momentum to the crust. 



The long history of glitches in Vela makes it possible to deduce some of 
the properties of the interior angular momentum reservoir independent of the 
details of the instability that triggers these events. Here we discuss the time 
distribution and average angular momentum transfer rate of Vela's ghtches and 
present evidence that ghtches in Vela represent a self -regulating instabihty for 
which the star prepares over a waiting interval. We obtain a lower limit on 
the fraction of the star's liquid interior responsible for glitches and discuss 
how this result can be used to constrain the dense matter equation of state and 
the structural properties of neutron stars. We conclude with discussion of the 
nearby isolated neutron star RX Jl 85635-3754, a promising candidate for a 
robust radius measurement that offers to constrain our understanding of dense 
matter and the origin of glitches. 

REGULARITY OF ANGULAR MOMENTUM TRANSFER 

A glitch of magnitude AOj requires angular momentum 



(1) 
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where Ic is the moment of inertia of the sohd crust plus any portions of the star 
tightly coupled to it. Most of the core liquid is expected to couple tightly to 
the star's solid component, so that Ic makes up at least 90% of the star's total 
moment of inertia [4]. [Glitches are driven by the portion of the liquid interior 
that is differentially rotating with respect to the crust]. The cumulative angular 
momentum imparted to the crust over time is 

J(t) = Icf^E^' (2) 

i 

where Vt = 70.4 rad s~^ is the average spin rate of the crust over the period of 
observations. Fig. 3 shows the cumulative dimensionless angular momentum, 
J(t ) / Ic^, over ~ 30 years of glitch observations of the Vela pulsar, with a linear 
least-squares fit. The average rate of angular momentum transfer associated 
with glitches is IcCtA, where A is the slope of the straight hne in Fig. 3: 

A = (6.44 ± 0.19) X 10~^ yr-\ (3) 

This rate A is often referred to as the pulsar activity parameter (see, e.g., [5]). 
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Figure 3 Cumulative dimen- 
sionless angular momentum, 
J/Ic^, imparted to the Vela 
pulsar's crust as a function of 
time. The straight line is a 
least-squares fit. 



The angular momentum flow is remarkably regular; none of Vela's 13 glitches 
caused the cumulative angular momentum curve to deviate from the hnear fit 
by more than 12%. Moreover, glitches occur at fairly regular time intervals; the 
standard deviation in observed glitch intervals is 0.53(Ai), where {At) = 840 
d is the average ghtch time interval. The probability of 13 randomly-spaced 
(Poisson) events having less than the observed standard deviation is only ~ 1%. 
These data indicate that Vela's glitches are not random, but represent a self- 
regulating process which gives a relatively constant flow of angular momentum 
to the crust with glitches occurring at fairly regular time intervals. 
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THE GLITCH RESERVOIR'S MOMENT OF INERTIA 

The frequent occurrence of large glitches requires that some fraction of the 
interior liquid (the reservoir) spins at a higher rate than the crust of the star. 
The average rate of angular momentum transfer in Vela's glitches (eq. 3) can 
be used to constrain the relative moment of inertia of the reservoir. Between 
glitches, the reservoir acquires excess angular momentum as the rest of the star 
slows under the magnetic braking torque acting on the crust. Excess angular 
momentum accumulates at the maximum rate if the reservoir coasts between 
glitches, without spinning down (Fig. 2). Hence, the rate at which the reservoir 
accumulates angular momentum capable of driving glitches is limited by 



«^res ^ -^resl^l) (4) 



where tl is the average spin-down rate of the crust, and /res is the moment of 
inertia of the angular momentum reservoir (not necessarily one region of the 
star). Equating Xes to the average rate of angular momentum transfer to the 
crust, IcO,A, gives the constraint. 



J^>-^A = G, (5) 

-'C i ' 



where the coupling parameter G is the minimum fraction of the star's moment 
of inertia that stores angular momentum and imparts it to the crust in glitches. 
Using the observed value of Vela's activity parameter A and = 22.6 Kyr, 
we obtain the constraint 



^ > Gvela = 1-4%. (6) 



A similar analysis for six other pulsars yields the results shown in Fig. 4. After 
Vela, the most significant limit is obtained from PSR 1737-30 which gives 
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Figure 4 The coupling pa- 
rameter G. The strongest con- 
straints are obtained for Vela 
and PSR 1737-30, for which 
13 and 9 glitches have been 
observed, respectively. Dia- 
monds indicate objects with 
only two observed glitches, 
for which error bars could 
not be obtained. References: 
0525-^21 [7], Crab [8], Vela 
[2], 1338-62 [9], 1737-30 [10, 
11], 1823-13 [11]. 

The similarity of G for tlie five objects of intermediate age suggests a common 
physical origin for glitches in these stars. The Crab pulsar and PSR 0525-1-21, 
however, appear to be unusual. It may be that the Crab's interior cannot ac- 
cumulate significant excess angular momentum between glitches, perhaps as a 
consequence of rapid thermal creep of superfluid vortices (see, e.g., [6]). The 
value of G for PSR 0525-1-21 is not well determined, since only two glitches 
from this object have been measured. 

IMPLICATIONS FOR THE DENSE MATTER EQUATION 
OF STATE 

The constraint of /res/^c ^ 1.4% for Vela applies regardless of where in the 
star glitches originate. A natural candidate for the angular momentum reservoir 
is the superfluid that coexists with the inner crust lattice [3], where superfluid 
vortex lines could pin to the nuclei and sustain a velocity difference between 
the superfluid and the crust. Within this interpretation, the constraint on /res/ Ic 
can be used to constrain the properties of matter at supranuclear density as we 
now demonstrate. 

The fraction of the star's moment of inertia contained in the soUd crust (and 
the neutron liquid that coexists with it) is given approximately by [12]: 

A/ _ 287r PtR"^ 
~T ~ 3 GM2 

Here nt is the density at the core-crust boundary, Pt is the pressure there, M and 
R are the stellar mass and radius, A = (1 — 2GM / Rc^)~^ is the gravitational 
redshift and m„ is the neutron mass. AI/I is a function of M and R with 
an additional dependence upon the equation of state (EOS) arising through the 
values of Pt and nt. Pt is the main EOS parameter as nt enters chiefly via a 
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correction term. In general, Pt varies over the range 0.25 < Pf < 0.65 MeV 
fm~^ for realistic equations of state [13]. Larger values of Pt give larger values 
for A///. 




Figure 5 Limits on the Vela pulsar's radius. The heavy dashed curve delimits allowed masses 
and radii that are compatible with the glitch constraint A//(/ — A/) > 1.4% for Pt — 0.65 
MeV fm""^ . The horizontal dashed lines indicate the mass limits for the survey of 26 radio pulsars 
of Ref. 16. For comparison, four mass-radius relations are shown for the equations of state of: 
Glendenning & Schaffner-Bielich (kaons, curve a) [14], Wiringa, Fiks & Fabrocini (curve b) 
[15], Glendenning & Schaffner-Bielich (kaons, curve c) [14] and a representative strange quark 
star (curve d). Thin curves are contours of constant radiation radius Roo- 

Combining eq. [7] with a lower limit on A/ and an upper limit on Pt gives 
a lower limit on the neutron star radius for a given mass. In order to relate our 
observational bound on Ircs/Ic to A/, we assume that the angular momentum 
reservoir is confined to the neutron superfluid that coexists with the nuclei of 
the inner crust. In this case, /res < and Ic > I — A/. Our observational 
limit on I^cs then gives AI/{I - A/) > AI/Ic > Ircs/h > 0.014. To obtain 
a conservative lower limit on the neutron star radius, we take Pt = 0.65 MeV 
fm~^ and nt = 0.075 fm~^. Eq. [7], gives the heavy dashed curve in Fig. 5. 
This curve is given approximately by 

i? = 3.6 + 3.9M/M0 . (8) 

Stellar models that are compatible with the lower bound on Ires must fall below 
this Une. For models above this line, the inner crust Uquid constitutes less than 
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1.4% of the star's moment of inertia, incompatible with the observed angular 
momentum requirments of the Vela pulsar. 

DISCUSSION 

Vela's mass and radius are unknown, so we cannot say where it falls in Fig. 5. 
However, mass measurements of radio pulsars in binary systems and of neutron 
star companions of radio pulsars give neutron star masses that are consistent 
with a remarkably narrow distribution, M = 1.35 ± 0.04Mq [16], indicated 
by the pair of horizontal dotted lines in Fig. 5. If Vela's mass is not unusual, 
eq. [8] constrains i?^ 8.9 km, under the assumption that ghtches arise in the 
inner crust superfluid. However, the quantity constrained by observations of 
the stellar luminosity and spectrum is not R but the larger 'radiation radius' 
Roo = A^/^R = (1 - 2GM/Rc^)-'^/'^R. If M = l.SBM© for Vela, the 
above constraint gives i?oo~ 12 km if glitches arise in the inner crust. For 
comparison, we show in Fig. 5 the mass-radius curves for several representative 
equations of state (heavy solid lines). Measurement of i?oo i$ 13 km would be 
inconsistent with most equations of state if M ~ l.SSM©. Stronger constraints 
could be obtained if improved calculations of nuclear matter properties give Pt 
significantly less than 0.65 MeV fm~^.^ For example, for M ~ I.SSMq, 
i^oo ~ 13 km would be required if Pt = 0.25 MeV fm~^. A measurement of 
i?oo ~ 11 km would rule out most equations of state regardless of mass or the 
angular momentum requirements of glitches, and could indicate that neutron 
stars are not made of neutrons at all, but of strange quark matter. Explaining 
ghtches in this case would be problematic, as strange stars have very thin crusts 
[17]. 

A black body fit to the unpulsed component of Vela's thermal emission gives 
Roo = 3 — 4 km [18]. This result is difficult to interpret without knowledge 
of the star's atmospheric composition and magnetic field strength; atmospheric 
effects could increase this estimate by a factor of up ~ 6, but could also decrease 
it by a factor ~ 2 [19]. Nevertheless, it would be interesting to check the extent 
to which our constraints on Vela are obeyed by other neutron stars. A promising 
candidate for a decisive measurement of a neutron star's radiation radius is RX 
J185635-3754, an isolated, non-pulsing neutron star [20]. A black body fit to 
the X-ray spectrum gives i?oo = 7.3(L'/120 pc) km where D is the distance 
(known to be less than 120 pc). Taken at face value, this result would not only 
be inconsistent with the radius requirements of ghtches in Vela - it would rule 
out all equations of state that do not involve strange matter. [It is possible that 



'if Pt is smaller than 0.65 MeV fm , the crust's moment of inertia would be smaller and the radius constraint 
more restrictive. For example, Pt = 0.25 MeV fm~^ moves the constraining contour to approximately 
R = i.7 + A.1M/Mq. 
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glitching pulsars are normal neutron stars, while RX Jl 85635-3754 is 'strange'] . 
However, either a non-uniform surface temperature or radiative transfer effects 
in the stellar atmosphere could raise this estimate significantly [21]. Recent 
HST observations of this source by F. Walter should give the proper motion 
and parallax, and hence, the distance. Future CHANDRA observations should 
yield more detailed spectral data and could establish the composition of the 
atmosphere if absorption lines are identified. If lines are present, atmospheric 
fits will give Roo, R and M, thus restricting RX Jl 85635-3754 to a region of 
Fig. 5. These data will undoubtedly further our understanding of matter at 
supranuclear density, and could establish whether neutron stars have properties 
consistent with an inner crust explanation of glitches. 

Acknowledgments 

We thank A. G. Lyne for providing us with glitch data for the Crab pulsar and P. M. McCuUoch 
for Vela data. This work was performed under the auspices of the U.S. Department of Energy, 
and was supported in part by NASA EPSCoR Grant #291748, NASA ATP Grants # NAG 53688 
and # NAG 52863, by the USDOE grant DOE/DE-FG02-87ER-40317, and by IGPP at LANL. 

References 

[1] A. G. Lyne, F G. Smith & R. S. Pritchard, Nature, 359, 
706 (1992). 

[2] J. M. Cordes, G. S. Downs & J. Krause-Polstorff, ApJ, 
330, 847 (1988); R M. McCulloch, et al. Aust. J. Phys., 
40, 725 (1987); C. Flanagan, lAU Circ. No. 4695 (1989); 
C. Flanagan, lAU Circ. No. 5311 (1991). 

[3] R W. Anderson & N. Itoh, Nature, 256, 25 (1975); M. 
Ruderman, 203, 213 (1976); D. Pines & M. A. Alpar, 
Nature, 316, 27 (1985). 

[4] M. A. Alpar, S. A. Langer & J. A. Sauls, ApJ, 282, 533 
(1984); M. Abney, R. I. Epstein & A. Olinto, ApJ, 466, 
L91 (1996). 

[5] A. G. Lyne, Lives of the Neutron Stars, p. 167. Ed: M. A. 
Alpar (Kluwer, 1995). 

[6] M. A. Alpar, P W. Anderson & D. Pines, ApJ, 276, 325 
(1984); B. Link, R. L Epstein, G. Baym, ApJ, 403, 285 
(1993); H. F Chau & K. S. Cheng, K. S., Phys. Rev. B, 

47, 2707 (1993). 

[7] G. S. Downs, ApJ, 257, L67 (1982). 

[8] P E. Boynton et al. ApJ, 175, 217 (1972); E. Lohsen, 
Nature, 258, 688 (1975); A. G. Lyne & R. S. Pritchard, 



Probing the Neutron Star Interior 9 

MNRAS, 229, 223 (1987); A. G. Lyne, R. S. Pritchard & 
F. G. Smith, Nature, 359, 706 (1992). 

[9] V. M. Kaspi, R. N. Manchester, S. Johnston, A. G. Lyne 
& N. D'Amico, ApJ, 399, L155 (1992). 

[10] A. G. McKenna & A. G. Lyne, Nature, 343, 349 (1990). 
[11] S. L. Shemar & A. G. Lyne, MNRAS, 282, 677 (1996). 
[12] B. Linli, R. I. Epstein & J. M. Lattimer, Phys. Rev. Lett., 
83, 3345 (1999). 

[13] J. M. Lattimer & M. Prakash, in preparation. 

[14] N. K. Glendenning & J. Schaffner-Bielich, Phys. Rev. C, 

60, 25803 (1999). 

[15] R. B. Wiringa, V. Fiks & A. Fabrocini, Phys. Rev., C, 38, 
1010 (1988). 

[16] S. E. Thorsett & D. Chakrabarty, ApJ, 512, 288. 

[17] C. Alcock, E. Farhi & A. Olinto, ApJ, 310, 261 (1986). 

[18] H. Ogelman, J. P. Finley & H. U. Zimmermann, Nature, 
361, 136 (1993). 

[19] M. Rajagopal, R. W. Romani & M. C. Miller, ApJ, 479, 
347 (1997). 

[20] F Walter, S. Wolk & R. Neuhauser, Nature, 379, 233 
(1996). 

[21] P. An, J. M. Lattimer, M. Prakash, & Walter, F. in prepa- 
ration. 



